Superfluidity of £ hyperons in /3-stable neutron star matter 
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In this work we evaluate the So energy gap of E~ hyperons in /3-stable neutron star matter. We 
solve the BCS gap equation for an effective E~E~ pairing interaction derived from the most recent 
parametrization of the hyperon-hyperon interaction constructed by the Nijmegen group. We find 
that the E~ hyperons are in a So superfluid state in the density region ~ 0.27 — 0.7 fm -3 , with a 
maximum energy gap of order 8 MeV at a total baryon number density of ~ 0.37 fm -3 and a E~ 
fraction of about 8%. We examine the implications on neutron star cooling. 
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Since the suggestion of Migdal pj, superfluidity in 
nuclear matter has received a great deal of attention 
over the last 40 years, partly due to its important 
consequences for a number of neutron star phenom- 
ena, such as pulsar glitches 0- 13- fl 01 and cooling 
rates 0, IH IS El- Nevertheless, whereas the pres- 
ence of superfluid neutrons in the inner crust of neutron 
stars, and superfluid neutrons together with supercon- 
ducting protons in their quantum fluid interior is well 
established and has been the subject of many studies 
[ISIIlllllliQEEOEIliita quantitative esti- 
mation of the pairing of other baryon species has not re- 
ceived so much attention up to date. In particular hyper- 
ons, which are expected to appear in neutron star matter 
at baryon number densities of order ~ 2no (no = 0.17 
fm~ 3 ), may also form superfluids if their interactions are 
attractive enough. It has been suggested that some neu- 
tron stars are cooled much faster than expected by a stan- 
dard cooling mechanism (i.e., modified URCA processes), 
and that more rapid and efficient mechanisms are needed 
H El El El Processes of the type Y B + 1 + n 
(e.g., A — > p + e~ + P e , £~ — > A + e~ + v e , etc) can pro- 
vide some of such rapid cooling mechanisms. Therefore, 
the study of hyperon superfluidity becomes of particu- 
lar interest since it could play a key role in them. The 
case of A superfluidity has been investigated by Balberg 
and Barnea (2J| using parametrized effective AA interac- 
tions. Results for A and S~ pairing using several bare 
hyperon-hyperon interaction models have been recently 
presented by Takatsuka et al. HE EE ■ The results of 
both groups indicate the presence of a A superfluid for 
baryon number densities in the range 2 — 4no. The latter 
authors suggest that both A and £~ become superfluid 
as soon as they appear in neutron star matter and that 
the formation of a £~ superfluid may be more likely than 
that of a A superfluid. 

Since the hyperon fraction (ny/ni,) in neutron star 
matter is not large (10% — 30% at most depending on the 
model), the Fermi momenta of hyperons are rather low 
although they appear at high values of the total baryon 



number densities. Therefore, the pairing interaction re- 
sponsible for hyperon superfluidity, if it exists, should 
be that due to the 1 5o wave which is most attractive 
at low momenta. In this work we evaluate the 1 S'o gap 
energies of £~ hyperons in /3-stable neutron star matter 
by solving the well-known BCS gap equation for an ef- 
fective pairing interaction derived from the most recent 
parametrization of the free baryon-baryon potentials for 
the complete baryon octet as defined by Stoks and Rijken 
[2^ |. We employ the model NSC97e of this parametriza- 
tion, since this model, together with the model NSC97f, 
results in the best predictions for hypernuclear observ- 
ables @. 

The crucial quantity in determining the onset of su- 
perfluidity is the energy gap function Afe, The value of 
this function at the Fermi surface is proportional to the 
critical temperature of the superfluid, and by determin- 
ing it we therefore map out the region of the density- 
temperature plane where the superfluid may exist. To 
evaluate it we follow the scheme developed by Baldo et 
al. p~E|. These authors introduced an effective pairing 
interaction Vk,w defined according to 
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which sums up all two-particle excitations above a cutoff 
momentum /cm > &f (&m = 2 fm -1 in this work). Previ- 
ous applications of this method to the neutron and pro- 
ton pairing [lfl have shown that it is stable with re- 
spect to variations of few, as we have al so confirmed. The 
quasiparticle energy E k is given by \J (e(k) — n) 2 + A^, 
being e(k) the single-particle energy in the medium for 
the particle species in question, [i the corresponding 
chemical potential, and Vk,k' the free baryon-baryon po- 
tential in momentum space, in our case the bare E _ S~ 
interaction of the NSC97e baryon-baryon potential. We 
note that the S~S~ channel is purely isospin 2 and there- 
fore there is no coupling to other hyperon-hyperon chan- 
nels in Eq. (QJ. For the 1 Sq channel the gap function can 



Baryon number density n b [fm ] 

FIG. 1: C omp osition of /3-stable neutron star matter. Taken 
from Ref. 

be determined by solving 

A fc = -£^'^;- (2) 

k' <kM 

Equations and Q are solved self-consistently and 
represent a totally equivalent formulation of the BCS gap 
equation. With this procedure: (i) a well-behaved pair- 
ing interaction is obtained, since the repulsive core of the 
bare interaction is integrated out and (ii) double count- 
ing of two-particle correlations is avoided. Excitations to 
intermediate states above ku are included in V, whereas 
excitations to states below kit are included in the gap 
equation (2J. We note here that for k > kp the domi- 
nant contribution to the quasiparticle energy Ek comes 
from the term (e(fc) — fi) 2 . Therefore, we can neglect 
in Eq. <QJ for k > k M > k F . Thus Eq. Q is decoupled 
from Eq. (01, and we can solve the linear equation for 
Vk.k' by the matrix inversion method before proceeding 
to solve the gap equation by iteration (see Ref. 0] for 
details). 

The relevant S~ fraction (shown in Fig. QJ, single- 
particle energy and chemical potential necessary to eval- 
uate Eqs. (QJ and {0| are taken from the Brueckner- 
Hartree-Fock calculations described in Ref. [3Jj], where 
the NSC97e baryon-baryon interaction was employed to 
describe the single-particle properties, the composition 
and Equation of State of /3-stable neutron star matter, 
and the neutron star structure. Therefore, to our knowl- 
edge, the present work is the first one which employs 
consistently the same baryon-baryon interaction model 
to determine the single-particle properties, the composi- 
tion, the Equation of State, the neutron star structure 
and the S~ energy gap. 

Figure El shows the energy gap Af of the S~ hyper- 
ons in /3-stable neutron star matter at T = with the 
composition shown in Fig. as a function of the total 
baryon number density. Although, as can be seen in Fig. 
01 the A may appear at higher densities, the 1 S AA 
matrix elements of the Nijmegen interaction (NSC97a-f) 
are all weakly attractive, and therefore the energy gap 



Total baryon number density n b [fm ] 

FIG. 2: Density dependence of the E~ energy gap Af in 
/3-stable neutron star matter at T = 0. 



for A hyperons is expected to be zero at all densities, i.e., 
these particles will unlikely form a superfluid within our 
model. This is at variance with the results of Balberg and 
Barnea [U Nevertheless, as stated before, these authors 
employed an effective parametrized interaction based on 
a G-matrix calculation to drive the gap equation and 
therefore overestimated, as pointed out by Takatsuka et 
al. 25, 26, 27], the A energy gap mainly due to double 
counting effects. Our results for the S~ are comparable 
to those of Takatsuka et al. j^E |22| which were obtained 
with several OBE hyperon-hyperon potentials. As these 
authors we find that S~ hyperons are in a 1 Sq superfluid 
state as soon as they appear in matter and that the S~ 
superfluid exists up to densities ~ 4no with a critical tem- 
perature T c ~ 10 10 K (see Fig.QJ. We find a maximum 
energy gap of order 8 MeV at a total baryon number den- 
sity of ~ 0.37 fm -3 and a E~ fraction of about 8%. This 
gap is quite large in comparison with the neutron and 
proton ones since the SE (and in particular the 
interaction in the Nijmegen NSC97a-f models is strongly 
attractive j^. We want to emphasize, however, that this 
strong attraction is questionable. Although these mod- 
els reproduce certain observables of A-hypernuclei, their 
predictions seem to be at odds with most of the scarce 
experimental data. The AA interaction, as mentioned 
before, is weak compared to the values deduced experi- 
mentally |3j|, and all types of hyperons are too strongly 
bound in nuclear matter This is especially suspect 
in the case of S~, since phenomenology of S~ atoms 
and hypernuclei j^J] indicate a much weaker, if not re- 
pulsive, S nuclear potential (see Ref. f° r a detailed 
discussion). Therefore, our results should be taken with 
caution. 

In Fig. we show the temperature dependence of the 
energy gap A F of S~ for several values of the total 
baryon number density and the corresponding /3-stable 
fractions of the S~. The gap function at finite tempera- 
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FIG. 3: Temperature dependence of the S~ energy gap Aj? 
in /3-stable neutron star matter. The fraction of S~ hyperon, 
n E - /n b , is indicated in each curve. The corresponding weak- 
coupling approximation (WCA) estimations for the critical 
temperatures are also indicated by the circle (nt = 0.3 fm -3 ), 
square (nt = 0.4 fm -3 ), diamond (rib = 0.5 fm -3 ) and trian- 
gle {n b = 0.6 fm" 3 ). 
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FIG. 4: Critical temperature of the So £ - superfluid as a 
function of the E~ number density. The internal temperature 
of evolved normal neutron stars is around 10 s K. 

ture can be obtained by solving 

A fe = - V VW^tanhf-^-) , (3) 
fc fc .* 2E k , \2k B T) ' v ; 

k'<k M 

where fcs is the Boltzmann's constant. We use the same 
approach as for the T = case. Here we ignore the tem- 
perature dependence in Vk,k' since for the temperature 
range of interest, ksT w — 4 MeV, the quasiparticle 
energy Ek for k > ku is at least of order 100 MeV, and 
thus we can ignore thermal excitations to states above 
fcM- In addition we use a "frozen" approximation for the 
single-particle energy, chemical potential and fraction of 



the X", i.e., we use the corresponding quantities obtained 
in the T — case, which is a reasonable approximation 
according to Refs. Therefore, as in the T = 

case, we first solve Eq. and then, with the effective 
interaction Vk,h' we solve Eq. Q). In Fig. we also show 
the critical temperatures estimated from the well-known 
weak-coupling approximation (WCA) 

k B T c « 0.57A F (T = 0) , (4) 

which is a reasonable good approximation as can be seen 
from the figure. 

Finally, in Fig. Q] we show the region in the 
temperature-X"-density plane where the E" hyperon is 
expected to be superfluid. Since the values of the critical 
temperature are all well above the typical internal tem- 
perature of evolved normal neutron stars (Ti nt ~ 10 8 K), 
the E - is in a 1 So superfluid state for number densities 
ranging from 2.3 x 10~ 4 fm -3 up to ~ 0.15 fm -3 , which 
corresponds, according to the composition shown in Fig. 

to a total baryon number density ranging from the X~ 
onset density (0.27 fm -3 ) to - 0.7 fm~ 3 (see Fig. 0. 

These results have implications for neutron star cool- 
ing. Since at low densities E~ is the only hyperon species 
that is present in our model, the most important con- 
tribution to the neutrino cooling rate at such densities 
comes from the reaction X - — > n + e~ + v e . In our 
model the threshold density for this reaction to occur 
is at around 1.6nrj. The direct action of such a rapid 
cooling mechanism, however, leads to surface tempera- 
tures much lower than that observed. Nevertheless, if the 
X~'s are in a superfluid state with energy gaps similar 
to what we found here, a sizeable reduction of the order 
cxp(— Af /ksT) may be expected in the neutrino emis- 
sivity of this process. Such a reduction will suppress the 
cooling rate and it will amount for neutron star surface 
temperatures more compatible with observation. Nev- 
ertheless, we should point out that this process will be 
also suppressed by the 3 P2 neutron pairing. This pair- 
ing exists practically for all supernuclear densities 
and, although it is relatively small (~ 0.1 MeV), it will 
suppress this process throughout almost the entire life of 
the neutron star. We want to finish just mentioning that 
hyperon superfluidity may be also important for r-mode 
stability calculations, since it may modify the tempera- 
ture and density dependence of hyperon bulk viscosity 
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